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ABSTRACT 

The existence of large-scale dynamos in rigidly rotating turbulent convection without shear is studied using 
three-dimensional numerical simulations of penetrative rotating compressible convection. We demonstrate 
that rotating convection in a Cartesian domain can drive a large-scale dynamo even in the absence of shear. 
The large-scale field contains a significant fraction of the total field in the saturated state. The simulation 
results are compared with one-dimensional mean-field dynamo models where turbulent transport coefficients, 
as determined using the test field method, are used. The reason for the absence of large-scale dynamo action 
in earlier studies is shown to be due to the rotation being too slow: whereas the a-effect can change sign, 
its magnitude stays approximately constant as a function of rotation, and the turbulent diffusivity decreases 
monotonically with increasing rotation. Only when rotation is rapid enough a large-scale dynamo can be 
excited. The one-dimensional mean-field model with dynamo coeflicients from the test field results predicts 
reasonably well the dynamo excitation in the direct simulations. This result further validates the test field 
procedure and reinforces the interpretation that the observed dynamo is driven by a turbulent a-efFect. This 
result demonstrates the existence of an a-effect and an a^-dynamo with natural forcing. 
Subject headings: MHD - convection - turbulence - Sun: magnetic fields - stars: magnetic fields 



1. INTRODUCTION 

Convective instability drives turbulence in the outer layers 
of late-type stars, such as the Sun. The large-scale magnetic 
fields of these stars are thought to arise from the interaction 
of turbulent convection and the overall rotation of the object. 
In mean-field dynamo theory (e.g. Moffatt 1978; Parker 1979; 
Krause & Riidler 1980; Rudiger & Hollerbach 2004), this pro- 
cess relies on large-scale differential rotation (O-effect) pro- 
ducing toroidal field by shearing and the a-effect which re- 
generates the poloidal field. In simple situations, the a-effect 
is related to the kinetic helicity of the flow (e.g. Steenbeck 
& Krause 1969). Dynamos where differential rotation is im- 
portant, e.g. the solar dynamo, are therefore often called aQ.- 
dynamos. 

The aO. dynamo process was first invoked by Parker (1955) 
to explain solar magnetism. Mean-field models have been 
used extensively ever since to study various aspects of dy- 
namos. Although very useful in their own right, these mod- 
els often rely on ill-known parameterizations of turbulence, 
such as the a-effect and turbulent diffusivity. Only during 
recent years numerical simulations have reached a level of 
sophistication where large-scale ffQ-dynamos have been ob- 
tained self-consistently within the framework of local Carte- 
sian simulations. These simulations operate in highly simpli- 
fied situations and the turbulence is driven by external forcing 
(e.g. Brandenburg et al. 2001; Brandenburg & Kapyla 2007; 
Kapyla & Brandenburg 2009). In rapidly rotating stars dif- 
ferential rotation is likely to play only a minor role (e.g. Hall 
1991) and there the a-effect generates also the toroidal field. 
These systems are called c^-dynamos. Again, such solutions 
have been found from direct simulations of forced turbulence 
(Brandenburg 2001; Mitra et al. 2009b,c). 

On the other hand, numerical simulations of magnetocon- 
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vection have been around for at least two decades, but some- 
what surprisingly large-scale dynamos were not found until 
quite recently (Jones & Roberts 2000; Rotvig & Jones 2002; 
Browning et al. 2006; Brown et al. 2007; Kapyla et al. 2008, 
hereafter Paper I; Hughes & Proctor 2009). The first two stud- 
ies are related to the geodynamo and therefore convection is 
rotationally dominated. The next two use a global models 
with an imposed tachocline (Browning et al. 2006) or rapid 
rotation (Brown et al. 2007), respectively. The last two are lo- 
cal models with imposed shear flows, reminiscent of the shear 
dynamos reported from non-helical turbulence (Yousef et al. 
2008a,b; Brandenburg et al. 2008). The origin of large-scale 
fields in the non-helically forced simulations cannot be due 
to the mean-field a-effect. However, the incoherent c-shear 
dynamo (e.g. Vishniac & Brandenburg 1997; Proctor 2007) 
or the mean-field shear-current effect (Rogachevskii & Klee- 
orin 2003, 2004; Kleeorin & Rogachevskii 2008) could drive 
these dynamos. In the shearing local convection simulations 
the same two effects have been suggested mainly due to the 
fact that rotating convection alone has so far been found un- 
able to generate large-scale magnetic fields (e.g. Nordlund et 
al. 1992; Brandenburg et al. 1996; Cattaneo & Hughes 2006; 
Tobias et al. 2008), although it should produce a mean a- 
effect according to theory. This drawback has prompted spec- 
ulations that the a-effect in its mean-field incarnation simply 
does not work (Hughes & Proctor 2009). 

A conclusive proof of the existence of an a^-dynamo should 
be the demonstration of large-scale dynamo action in a sim- 
ulation without shear An additional property of such a setup 
is that the Vishniac & Cho (2001) flux of magnetic helicity is 
absent. This might lead to slow saturation of the large-scale 
magnetic field, unless there are other ways of shedding small- 
scale magnetic helicity. 

The most natural way of explaining the large-scale mag- 
netic fields seen in the rapidly rotating global convection sim- 



2 



ulations (Brown et al. 2007) would be in terms of a mean-field 
a^- or aQ-dynamo. In the global models large-scale shear 
flows can also be generated, so it is not straightforward to de- 
termine what type of dynamo is operating there. In the present 
study, we use local simulations to demonstrate that large-scale 
dynamos can indeed be excited in rapidly rotating convec- 
tion, i.e. in the absence of shear, as long as rotation is rapid 
enough. We determine the turbulent transport coefficients 
from the simulations using the test field procedure (Schrin- 
ner et al. 2005, 2007) and use them in a one-dimensional dy- 
namo model to test their consistency with the direct simula- 
tions. A more thorough study of the turbulent transport coeffi- 
cients from convection using the test field method is presented 
elsewhere (Kapyla et al. 2009, hereafter Paper 11). 

2. THE MODEL 

Our model setup is similar to that used by Brandenburg 
et al. (1996), Ossendrijver et al. (2001, 2002), Kapyla et al. 
(2004, 2006), and those in Papers 1 and 11. A rectangular por- 
tion of a star is modeled by a box situated at colatitude 6. 
The box is divided into three layers, an upper cooling layer, a 
convectively unstable layer, and a stable overshoot layer (see 
below). The following set of equations for compressible mag- 
netohydrodynamics is being solved: 
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where D/Dt - d/dt +U ■'V. The magnetic field is written in 
terms of the magnetic vector potential. A, with JB = V x A, 
J = ;Uq ' V X B is the current density, fio is the vacuum perme- 
ability, 77 and V are the magnetic diffusivity and kinematic vis- 
cosity, respectively, K is the heat conductivity, p is the density, 
U is the velocity, g - -gz is the gravitational acceleration, 
and ft = f2o(-sin0, 0, COS0) is the rotation vector. The fluid 
obeys an ideal gas law p - iy - l)pe, where p and e are pres- 
sure and internal energy, respectively, and y - cp/cy - 5/3 
is the ratio of specific heats at constant pressure and volume, 
respectively. The specific internal energy per unit mass is re- 
lated to the temperature via e - cyT. The rate of strain tensor 
S is given by 

Sij = ^(Uij + Ujj) - ■ U. (5) 

The last term of Eq. (|4|i describes cooling at the top of the do- 
main. Here t(z) is a cooling time which has a profile smoothly 
connecting the upper cooling layer and the convectively un- 
stable layer below, where t — > 00. 

The positions of the bottom of the box, bottom and top of 
the convectively unstable layer, and the top of the box, respec- 
tively, are given by (z\,Z2,Z3,Z4) = (-0.85,0, 1, 1.15)(i, where 
d is the depth of the convectively unstable layer. Initially the 
stratification is piecewise polytropic with polytropic indices 
{m\,m2,m^) - (3, 1, 1), which leads to a convectively unsta- 
ble layer above a stable layer at the bottom of the domain. The 
cooling term leads to a stably stratified isothermal layer at the 
top. The horizontal extent of the box Lh = - L, is varied 
between 2d and %d. All simulations with rotation use = 0° 
corresponding to the north pole. 



2. 1. Units and nondimensional parameters 
Nondimensional quantities are obtained by setting 

<^ = ^ = PO = Cp = yUo = 1 , 



(6) 



where po is the initial density at zi- The units of length, time, 
velocity, density, entropy, and magnetic field are 



[x]^d, [t]^y[drg, [f/] = v^, bo]=po 

[s] = Cp , [B] = ^JdgpopQ . 



(7) 



We define the fluid and magnetic Prandtl numbers and the 
Rayleigh number as 



V V 
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where ;t'o = ^/(PmCp) is the thermal diffusivity, and p^ is the 
density in the middle of the unstable layer, Zn, - jiz?, - 22)- 
The entropy gradient, measured at Zm, in the nonconvecting 
hydrostatic state, is given by 



]_ds\ ^ V-Vad 

cpdz/o Hp 



(9) 



where V - Vad is the superadiabatic temperature gradient 
with Vad = l-l/y, V - (dlnT /din p)~^, and where Hp 
is the pressure scale height (Brandenburg et al. 2005). The 
amount of stratification is determined by the parameter ^0 - 
(y - l)eQ/{gd), which is the pressure scale height at the top of 
the domain normalized by the depth of the unstable layer We 
use ^0 - 1 /3 in all cases, which results in a density contrast of 
about 23 across the domain. We define the fluid and magnetic 
Reynolds numbers via 

Re = ^, Rm=^=PmRe, (10) 
vk{ rjkf 

where kf - 2jT/d is adopted as an estimate for the wavenum- 
ber of the energy-carrying eddies. Note that with our defi- 
nition Rm is smaller than the usual one by a factor 2/:. The 
amount of rotation is quantified by 



Co = 



2Qo 

^rms^f 



(11) 



In order to facilitate comparison with earlier studies of rotat- 
ing Rayleigh-Benard convection we also quote the value of 
the Taylor number. 



Ta = (2Qot/Vv) = (47r2CoRm/Pm) . 
The equipartition magnetic field is defined by 

Beq = {popU^y^\ 

where angular brackets denote volume averaging. 



(12) 



(13) 



2.2. Boundary conditions 

Stress-free boundary conditions are used in the vertical (z) 
direction for the velocity, 

t/^_, = f/,,, = f/- = 0, (14) 

where commas denote partial derivatives, and either vertical 
field of perfect conductor conditions are used for the magnetic 
field, i.e. 

B,_ = B,, = (vertical field), (15) 
Bx,z-By^7-B--Q (perfect conductor). (16) 
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TABLE 1 

Summary of the runs. Here, Ma = Uj^s/Cgi^)''^, firms = firms/Beq: and the B'*' are the sum of the rms-values of the Fourier amplitudes of the horizontal 

COMPONENTS OF THE MAGNETIC FIELD FOR MODE k NORMALIZED BY Bnils- VALUES WITH S'"' > 0.07 OR B*" > 0.15 ARE SHOWN IN BOLD FACE AND DEDICATE THE PRESENCE OF 
SIGNIFICANT LARGE-SCALE FIELDS. ThE LAST COLUMN DENOTES THE MAGNETIC FIELD BOUNDARY CONDITION AT THE Z-BOUNDARIES. 



Run 


grid 




Pr 


Pm 


Ra 


Rm 


Co 


Ta 


Ma 


- 

firms 


5(0) 


5(1) 


bC 


Al 


128^ 


2 


0.27 


2 


1.5 ■ 


10*^ 


70 








0.044 


0.12 


0.02 


0.05 


vf 


A2 


128' 


2 


0.27 


2 


1.5 ■ 


10*^ 


68 


0.37 


2.5 ■ 10^ 


0.043 


0.23 


0.04 


0.10 


vf 


A3 


128' 


2 


0.27 


2 


1.5 ■ 


lO*" 


68 


0.74 


10* 


0.043 


0.25 


0.04 


0.09 


vf 


A4 


128' 


2 


0.27 


2 


1.5 ■ 


lO*" 


69 


1.5 


4- 10* 


0.043 


0.14 


0.02 


0.05 


vf 


A5 


128' 


2 


0.27 


2 


1.5 ■ 


10* 


62 


4.1 


2.5 ■ 10' 


0.039 


0.18 


0.04 


0.07 


vf 


A6 


128^ 


2 


0.27 


2 


1.5 ■ 


10** 


42 


11.6 


108 


0.028 


0.38 


0.22 


0.18 


vf 


A7 


128^ 


2 


0.27 


1 


1.5 ■ 


lO*" 


22 


11.6 


108 


0.028 


0.27 


0.15 


0.16 


vf 


A8 


128^ 


2 


0.27 


0.67 


1.5 ■ 


lO*" 


17 


9.8 


lO^ 


0.033 


0.00 


0.00 


0.00 


vf 


A9 


128^ 


2 


0.27 


0.4 


1.5 • 


10** 


10 


9.7 


108 


0.033 


0.00 


0.00 


0.00 


vf 


AlO 


2563 


2 


0.14 


2 


3.0- 


10* 


107 


9.5 


4- 10* 


0.034 


0.41 


0.08 


0.15 


V/ 


Bl 


256^ X 128 


4 


0.27 


2 


1.5 ■ 


10* 


77 








0.048 


0.22 


0.02 


0.06 


vf 


B2 


256^ X 128 


4 


0.27 


2 


1.5 ■ 


10* 


67 


0.38 


2.5 • 10^ 


0.042 


0.25 


0.02 


0.05 


vf 


B3 


256^ X 128 


4 


0.27 


2 


1.5 • 


10* 


66 


0.76 


10* 


0.042 


0.22 


0.02 


0.05 


vf 


B4 


256^ X 128 


4 


0.27 


2 


1.5 • 


10* 


67 


1.5 


4- 10* 


0.042 


0.16 


0.01 


0.03 


vf 


B5 


256^ X 128 


4 


0.27 


2 


1.5 ■ 


10* 


55 


4.6 


2.5 ■ 10' 


0.035 


0.31 


0.05 


0.09 


vf 


B6 


256^ X 128 


4 


0.27 


2 


1.5 • 


10* 


42 


11.9 


lO* 


0.027 


0.41 


0.12 


0.15 


vf 


B7 


256^ X 128 


4 


0.27 


2 


1.5 • 


10* 


41 


12.4 


lO** 


0.026 


0.48 


0.18 


0.26 


pc 


CI 


512-^ X 128 


8 


0.27 


2 


1.5 ■ 


10* 


42 


12.0 


10** 


0.026 


0.45 


0.07 


0.17 


vf 


Dl 


512^x256 


4 


0.14 


2 


3.0- 


10* 


97 


10.5 


4 ■ 10** 


0.030 


0.58 


0.08 


0.30 


vf 



The vertical field (VF) conditions permit a magnetic helicity 
flux, but no Poynting flux, whereas the perfect conductor (PC) 
conditions do not allow helicity fluxes either. In the x and y 
directions periodic boundary conditions are used. The sim- 
ulations were made with the Pencil Code', which uses sixth 
order explicit finite diff'erences in space and third order accu- 
rate time stepping method. Resolutions of up to 512^ x 256 
mesh points were used. 

3. RESULTS 

We perform a parameter study where we vary the system 
size, the effect of rotation, quantified by Co, and the Reynolds 
numbers in order to study the existence of large-scale dy- 
namos in rotating convection. The runs are listed in Table [1] 
We briefly describe the hydrodynamics of these runs in § 13.11 
and the results on dynamo excitation and large-scale magnetic 
fields are given in § 13.21 Finally, in § l3.3l and § l3.4l the turbu- 
lent transport coefficients and the interpretation of the results 
in the framework of mean-field electrodynamics, respectively, 
are presented. 

3.1. The hydrodynamic state 

All our simulations start with a small (of the order of 
lO^^Beq) random magnetic field. This field remains dynami- 
cally insignificant in the initial stages, which, in most cases, 
span the first 200 turnover times of the run. As in Paper I we 
consider this interval to represent the hydrodynamic state of 
the simulation. 

The effect of rotation on convection has been described in 
various papers (e.g. Kapyla et al. 2004; Giesecke et al. 2005, 
Cattaneo & Hughes 2006; Hughes & Cattaneo 2008; Tobias 
et al. 2008). The most visible effect is the decreasing size of 
convective cells as rotation is increased; see Fig. [T] Without 
rotation (Run Bl) and for horizontal extent Lu/d = 4 the flow 
is dominated by essentially a single large cell, but already for 
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Co X! 1 .5 (Run B4) there are of the order of ten or more cells 
in the domain. When rotation is increased further, the cell 
size continues to diminish. This feature is familiar from sim- 
ulations of unstratified convection (King et al. 2009), although 
there the flow appears to become more easily laminar at large 
Taylor numbers (Ta = 10**) and modestly large Rayleigh num- 
bers (Ra = 1.5 • 10^) that we used here (see Table[T]i. 

The decrease in cell size is also manifested by two- 
dimensional power spectra of the velocity, see the upper panel 
of Fig. |2] For the nonrotating case (Run Bl) the most power 
is found at k/ki = 1, where = In/Ln. There is a tendency 
for the wavenumber of the maximum, kmux, to shift towards 
higher k for more rapid rotation. The lower panel of Fig. |2] 
shows the velocity power spectra for Runs A6, B6, and CI 
with Lii/d = 2, 4, 8, respectively, shifted so that the same spa- 
tial scales coincide. All runs have Re ^ 28 (21) and Co ^ 9 
(12) in the kinematic (saturated) regime. For the smallest sys- 
tem size the energy still peaks near the box scale. When the 
box size increases the spatial scale at which most energy is 
found stays the same and, as a consequence, the scale separa- 
tion, ^max/^i, between the energy carrying scale and the box 
scale increases. 

Our finding that A:max is independent of the horizontal sys- 
tem size is important because large-scale dynamos require 
some amount of scale separation, i.e. that the power spec- 
trum of the turbulence should peak at some scale smaller than 
the system size. This allows the energy to cascade to smaller 
and larger scales. The former eventually dissipates the energy 
into heat whereas the latter can generate large-scale structures 
through an inverse cascade of energy. The turbulent a-effect 
can also be interpreted as an inverse cascade (Brandenburg 
2001). If, however, the flow energy peaks at, or close to, the 
scale of the system, it is impossible to cascade energy to larger 
scales. 



3.2. Dynamo excitation and large-scale magnetic fields 
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Fig. 1. — Vertical velocity U- from the middle of the convectively unstable layer from Runs B1-B6 in the kinematic stage at ? = 400 -\jd/g, corresponding to 
tu,nisk[ » 70 ... 1 10, depending on the run. Light (dark) color indicates ascending (descending) motion. 



Recent numerical studies have indicated that shear plays an 
important role in exciting large-scale dynamos (e.g. Yousef et 
al. 2008a,b; Paper I; Hughes & Proctor 2009). Although the 
process generating the large-scale magnetic field in the non- 
helical case is still under debate (e.g. Rogachevskii & Klee- 
orin 2003, 2004; Rudiger & Kitchatinov 2006; Proctor 2007; 
Brandenburg et al. 2008; Schekochihin et al. 2009; Sridhar & 
Subramanian 2009), there are some indications that in helical 
turbulence with shear a classical aO., or ff-shear dynamo, si- 
multaneously with a shear-current and ft x J-dynamos, might 
be the explanation (Paper II). However, in the absence of 
shear, large-scale dynamos due to helical turbulence have so 
far been obtained only in idealized systems where the helical 
flow is driven by external forcing (Brandenburg 2001; Mi- 
tra et al. 2009b, c). In particular, the lack of large-scale dy- 
namos in rotating convection has been puzzling (e.g. Cattaneo 
& Hughes 2006), although an c-effect should be present ac- 
cording to theory. One possible explanation is that, whereas 
the turbulence in more idealized studies can be almost fully 
helical owing to the forcing, in rotating convection the frac- 
tional helicity is often much smaller (e.g. Brandenburg et al. 
1990, 1996; Kapyla et al. 2004). 

In Paper II we found that the turbulent a-effect increases 
and turbulent diffusivity decreases, respectively, as functions 
of rotation for small fluid Reynolds numbers (see below). 
These results suggest that a large-scale dynamo in rotating 
convection is excited only if rotation is rapid enough. Our 



conjecture is that this regime might not have been reached in 
earlier studies thus failing to produce large-scale dynamos. 

The evolution of the rms-value of the total magnetic field 
for Runs B1-B6 is shown in Fig. [3] We find that the growth 
rate of the field increases for Co < 0.74 (Runs B1-B3) after 
which it starts to decrease. However, for the most rapid rota- 
tion (Run B6) the growth rate again increases and is close to 
the nonrotating case. The saturation level of the magnetic field 
is practically constant for Co < 0.74 (Runs B 1-B3) and some- 
what lower for Co a; 1.5 (Run B4). For the two most rapidly 
rotating cases the saturation level is generally higher and more 
variable. In these runs the maxima of the field are associated 
with periods where significant amounts of large-scale fields 
are present; see Fig. H] for the time evolution of horizontally 
averaged fields By and By for Run A6. Comparing the field 
in the kinematic and saturated stages (Fig. |5]l shows that in 
the kinematic regime the field is concentrated in small scales, 
whereas in the late stages a clear large-scale structure is vis- 
ible. The behavior of the growth rate and saturation level of 
the field as functions of rotation could be understood as fol- 
lows: the small-scale dynamo is enhanced for slow rotation 
but starts to be rotationally quenched for Co > 1 . This might 
be explicable by an associated decrease in the length scale of 
the turbulence. On the other hand, for rapid enough rotation, 
i.e. Co > 3, the large-scale dynamo becomes excited and in- 
creases the growth rate and saturation level. We stress that the 
proposed suppression of the small-scale dynamo action due 





1 10 100 

fc/fcl 



Fig. 2. — Two-dimensional power spectra of velocity from the middle of the 
convectively unstable layer as functions of rotation (upper panel) and system 
size (lower panel) in the kinematic regime at f = 400 ^jd|g. Power law slopes 
proportional to kT^^^ are shown for reference. The spectra in the lower panel 
have been shifted so that the same spatial scales coincide on the .v-axis and 
scaled by Ln/d so that the curves lie on top of each other 
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Fig. 3. — Upper panel: root mean square of the total magnetic field as a 
function of time for the same runs as in the upper panel of Fig. [2] Lower 
panel: the same as above but in linear scale. 



to rotation is at this stage only a conjecture that is consistent 
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Fig. 4. — Horizontally averaged magnetic field components (top 
panel) and By (bottom panel), normalized by the equipartition magnetic field 
strength. From Run A6 with Co « 12 and Rm x 42. The horizontal dotted 
lines at z/d = and z/d = 1 indicate the base and top of the convectively 
unstable layer, respectively. 



with the numerical data. Whether this interpretation is correct 
should be studied in more detail separately. It should also be 
noted that we cannot expect rotation to quench the small-scale 
dynamo in the astrophysically relevant regime where Rm is 
much larger than in the present simulations. There are, how- 
ever, no compelling reasons to expect that the large-scale dy- 
namo is suppressed when a small-scale dynamo is also present 
as long as magnetic helicity fluxes are allowed to leave the do- 
main. For the most rapidly rotating case, with Co 9 in the 
kinematic regime, the critical magnetic Reynolds number is 
roughly 20, see Fig. |6] In Paper I we reported that in the ab- 
sence of rotation, the critical Rm for the excitation of a small- 
scale dynamo is roughly 30 for a similar system. 

In Paper I we found the large-scale field to be dominated by 
the ki- = fc, = mode and it could thus be well described by 
a horizontal average. However, in the present case the mean 
fields show a more complicated structure and tend to have 
comparable contributions from A;^, ky + modes as well (see, 
e.g. the right panel of Fig.|5]l. Thus a simple horizontal aver- 
age no longer suflices to represent the large-scale field. For a 
better measure, we perform a two-dimensional Fourier trans- 
form of the field at each horizontal layer, the result of which 
we denote by Bi{kx,ky,z,t). Now we project this quantity 
onto a one-dimensional wavenumber k^ = k^ + k^ and denote 

its Fourier amplitude by B^-"^\z, t), where m is the discretized 
wavenumber bin corresponding to k/ki. So, for k/ki = 0, 
Bf \z, t) is the same as a horizontal average B,. For larger k, 

Bf^\z, t) gives a measure of the strength of the mode k. In 
the following we restrict the analysis to the smallest values of 
k/ki = ... 2 that describe the large-scale field. 

Representative results for the horizontal components of 
Bf \z, t) and Bf\z, f) from Run B6 with Co =s 12 are shown 
in Fig. |2] The large-scale field seems to show opposite signs 
in the convectively unstable and stable layers. Sign changes 
also occur but they are rather irregular and do not appear to 
follow a consistent cycle. The root mean square of the ampli- 
tudes for the kjki - and 1 contributions to the horizontal 
magnetic fields are of the order of ten per cent of the equipar- 
tition value. This is to be contrasted with Fig. [8] where the 
sums of the rms-values of the amplitudes of the three smallest 





10" 




Rm~56 (Run A8) 


10-= 




y Rm~28 (Run A7) 


10-" 


















10-° 




■■^v^s Rni~18 (Run A8) 






Rm~ll (Run A9) 



200 400 600 800 1000 



Fig. 6. — Root mean square total magnetic field as a function of time for 
four magnetic Reynolds numbers for Runs A6-A9. 



wavenumbers are shown as functions of rotation. Substantial 
large-scale fields are observed only for the two largest values 
of Co; for Co < 1.5 the runs show very weak large-scale con- 
tributions. This is also visible in the two-dimensional power 
spectra, taken from the middle of the convectively unstable 
layer; see Fig.|9] The most rapidly rotating case (Run B6) is 
the only one showing clear signs of large-scale fields, in ac- 
cordance with the fact that the large-scale field is only period- 
ically present in Run B5. The velocity spectra in the saturated 
state are similar to those in the kinematic phase. 

In comparison to earlier studies of rotating convection (e.g. 
Cattaneo & Hughes 2006; Tobias et al. 2008), we note that it is 
characteristic of these studies that when the Rayleigh number 
is increased, the Taylor number is kept constant. Increasing 
Ra in these models generates a larger Mrms and this inevitably 
means that the Coriolis number decreases as the Rayleigh 
number is increased, i.e. for a fixed Ta the rotational influence 
is large for a small Rayleigh number and vice versa. For ex- 
ample, in the paper of Cattaneo & Hughes (2006), the smallest 
Rayleigh numbers in combination with Ta = 5 ■ 10^ gives a 
Coriolis number comparable to our largest values. However, 
these simulations do not exhibit dynamo action due to a too 
low Rm, whose value also depends on M„ns. For their high- 
est Rayleigh number case, however, Rm is large enough for 
dynamo excitation but the Coriolis number is smaller by ap- 



proximately an order of magnitude and no large-scale fields 
are observed. Similar arguments apply to the simulations of 
Tobias et al. (2008). 

We find that the large-scale dynamo is excited for all box 
sizes for the most rapidly rotating case explored in the present 
study, as is evident from the spectra shown in Fig.[TO] From 
the spectra it would seem that an increasing amount of energy 
is in the large scales as the box size increases. The growth 
rate of the total field does not show any clear trend with the 
system size: the largest departure from a constant growth rate 
is the somewhat lower value for Run B6 with the intermediate 
box size (see Fig.fTTTi. 

Figure [TT] shows the sums of the Fourier amplitudes of the 
three smallest wavenumbers as functions of the system size. 
For the smallest box (Run A6), the large-scale field is more 
concentrated on the kjkx - contribution, whereas in Run B6 
with L-a_ld - 4 the amplitudes for kjk\ =0 and 1 are similar. 
For the largest domain size the = mode is significantly 
weaker than the kjk\ - 1 and 2 modes. These results suggest 
that for the present parameters the maximum size of the large- 
scale structures is somewhere in the range 2 < Lmax/t/ < 8. 

Comparing the saturation level of the large-scale magnetic 
field in the small box Runs A7, A6, and AlO shows a signif- 
icant decrease in the m = component in Run AlO whereas 
the strength of the m - 1 mode is only mildly affected. On 
the other hand, comparing Runs B6 and Dl with a larger do- 
main size shows again a decreasing m = contribution in 
Run Dl but a two times larger m - \ component. However, 
these numbers should be taken only as a rough guide because 
the large-scale contribution to the magnetic field shows large 
fluctuations and the higher Rm runs are fairly short. Taken at 
face value, the results would seem to suggest that the strength 
of the m - mode decreases with increasing Rm and that 
the m - I mode remains unaffected or that it can even in- 
crease. We note that the highest Rm runs also have larger 
fluid Reynolds and Rayleigh numbers which means that also 
the flow is more turbulent in those cases which could affect the 
dynamo and thus the saturation level of the large-scale field. 

Although we have used open (vertical field) boundary con- 
ditions that do permit magnetic helicity fluxes, such fluxes 
may not actually occur unless they are driven toward the 
boundaries by internal magnetic helicity fluxes. One such flux 
is the Vishniac-Cho (2001) flux, but it requires shear which 
is absent in our case. Other fluxes are possible (Subrama- 
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From top to bottom: Bf \z, t), b'."(z, t), their rms-values; Bf \z, t), b\^\z, t), 
and their rms-values. 



nian & Brandenburg 2006), but we do not know how efficient 
they are in our model. It is therefore unclear whether one 
should expect the saturation of the large-scale field to occur 
on a dynamical or a resistive time scale, and what the rele- 
vant length scale is. In Fig. [12] we show the saturation be- 
havior for Run B6 and compare with a curve proportional to 
1 - exp[-2rik^{t - fs)] that would be expected for resistively 
dominated saturation behavior (Brandenburg 2001). Here, 
k - kf has been chosen and ts marks the end of the linear 
growth phase, which is also the time when the small-scale 
magnetic field has saturated. The result is not entirely conclu- 
sive, and larger magnetic Reynolds number would be needed 
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Fig. 8. — The sums of the root mean square values of the Fourier amplitudes 
of B v and By for the modes k/k[ = ... 2 as function of rotation for the 
Runs B1-B6. Linestyles as in the upper panel of Fig.|2] 




Fig. 9. — Two-dimensional power spectra of velocity (upper panel) and 
magnetic field (lower panel) from Runs B1-B6 in the saturated state. Power 
laws proportional to k^^^^ are shown for reference. 



to clarify this, but it is certainly possible that saturation of the 
large-scale field is resistively dominated. Similar behavior is 
certainly expected in the case of perfect conductor boundary 
conditions (Run B7; see the dashed lines in Fig.[T2li. We note 
that the saturation level of the total and mean fields are higher 
in Run B7 with perfect conductor boundary conditions than in 
Run B6 with vertical field conditions. The behavior is quali- 
tatively similar to the forced turbulence simulations of Bran- 
denburg (2001). 

3.3. Turbulent transport coefficients 

We attempt to connect the dynamos seen in the direct simu- 
lations to theoretical considerations by comparing the results 
to a mean-field dynamo that relies on turbulent transport co- 
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Fig. 10. — Two-dimensional power spectra of velocity (upper panel) and 
magnetic field (lower panel) as functions of system size. Linestyles and scal- 
ing as in the lower panel of Fig. |2] In the lower panel the upper curves show 
the spectra from the saturated state whereas the lower curves show the spectra 
from the kinematic state multiplied by 10^. 

efficients, such as the a-effect and turbulent diffusivity. In 
order to extract these coefficients we use the test field method 
(Schrinner et al. 2005, 2007) in the kinematic regime. De- 
tailed descriptions of the method, as implemented here, have 
been presented elsewhere (e.g. Brandenburg et al. 2008; Sur 
et al. 2008; Mitra et al. 2009a; Paper II); here we only outline 
the general procedure. The method relies on a set of orthogo- 
nal test fields that do not evolve in time and do not react back 
onto the flow. These properties allow an unambiguous de- 
termination of the turbulent transport coefficients in the kine- 
matic regime without the complications associated with some 
other methods that can be used. 
We restrict the study to test fields that depend only on z 

B'' = Bo(cosfe,0,0), B^' = Bo(0,cosfe,0), (17) 
B^' = Bo(sinfe,0,0), B^' = Bo(0, sinfe, 0), (18) 

where Bq is the amplitude of the field and k is the wavenumber 
of the test field. (The results are strictly independent of the 
value of so in non-dimensional units it can be set to one.) 
We use kjki - 1, where - 2nlL,, in all models. The k- 
dependence of the coefficients in convection was studied in 
Paper II. The electromotive force can be written as 

£; = aijBj - }]ijfioJj, (19) 

where 77,1 - 77,23 and 77,2 = -77,13 • The 4+4 coefficients are 
then obtained by inverting a simple matrix equation, relating 
the rank-2 tensor components to rank-3 tensor components. 

It is convenient to discuss the results in terms of the quanti- 
ties 

q;= j(au + a'22), r = 5(«2i - an), (20) 
'7t= 5(^/11 + ?722), ^ = 5(?72i - '?12)- (21) 
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Fig. 11. — Top panel: rms- values of the total magnetic field as functions of 
the system size. The three lower panels show the sums of the rms-values of 
the Fourier amplitudes of B , and B,. for ^/^i = (second panel from the top), 
kjkx = 1 (second from the below), and k/ki = 2 (bottom panel). Linestyles 
as in the lower panel of Fig. [2] 



These quantities can be understood to represents the genera- 
tion of large-scale magnetic fields (a), turbulent pumping of 
magnetic fields (y), turbulent diffusivity (771), and the SI x J 
effect (6). To normalize our results, we use isotropic expres- 
sions derived under the first order smoothing approximation 
(see, e.g. Paper II) 

ffO = |"rms, ?7tO = ^Mrms^f (22) 

Figure [13] shows the kinetic helicity a; • u, where u) — 
V X It, and the turbulent transport coefficients for Runs B2- 
B6 as functions of rotation. The magnetic Reynolds number 
varies between Rm ^ 58 ... 86 with larger values occurring 
for slower rotation. Pm = 2 used in all runs. These simula- 
tions differ from those in Paper II in that Pr and Pm are respec- 
tively five times and 2.5 times smaller We find that the kinetic 
helicity increases monotonically as a function of rotation as 
in Paper II, although here the increase is not as steep as in the 
previous results. The (T-effect is approximately constant for 
Co X 0.3 ... 1.2, whereas for Co w 3.3 the positive region in 
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the top half of the convectively unstable layer disappears giv- 
ing nearly zero value. For the most rapid rotation, Co ^ 8.7, 
a is negative in the whole convection zone. At first sight this 
result seems to contradict the corresponding results of Paper 
II, but we find that there is a qualitative change in the behav- 
ior of a in the rapid rotation regime when the fluid Reynolds 
number is increased. For Re 2 there is no sign change and 
a monotonously increasing magnitude whereas for Re w 7 the 
positive region in the upper layers of the convectively unsta- 
ble layer disappears. Increasing the Reynolds number further 
to roughly 30 gives the result shown in Fig. [13] We note that 
controlling Re and Co a priori is diflicult in the rapid rotation 
regime because convection is increasingly suppressed espe- 
cially for small Reynolds numbers. 

For the slowest rotation, the turbulent pumping, y, is very 
similar in profile and magnitude as in Paper II. In the case of y 
the behavior seems again different with negative values in the 
whole convection zone for Co > 1.17. The magnitude of y is 
significantly decreased for our case with the most rapid rota- 
tion. The turbulent diffusivity is rotationally quenched simi- 
larly as in Paper II. Also the coefficient S behaves in a similar 
fashion as in Paper II. 

The relatively unchanged magnitude of a and the signifi- 
cantly reduced value of rii suggests that excitation of a mean- 
field dynamo could be possible in our cases with the most 
rapid rotation. One way to quantify this is to compute a local 
dynamo number 

caiz) = — , (23) 

where t/t = ?7t + 1- The results are shown in Fig. [14] We find 
that the dynamo number peaks in the overshoot layer just be- 
low z/d = because the profile of the a-effect extends some- 
what deeper than that of turbulent diffusivity. The magnitude 
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layer, respectively. 



of the maximum of Cq, in the overshoot increases monotoni- 
cally until Co w 3.3 but the value in the convectively unstable 
layer remains small. For the most rapidly rotating case, Ca(z) 
increases significantly also in the convectively unstable region 
due to the rotationally quenched value of t]^. In the homoge- 
neous case, where c„ is constant, dynamo action is possible 
when its value based on the lowest wavenumber in the domain 
exceeds unity, i.e. Ca = CakfL,/27T > 1. In our case, owing 
to the presence of the overshoot layers, kfL./ln: - 2, so one 
might expect that dynamo action is possible when the average 
value of Caiz) exceeds 0.5. For our case with the largest rota- 
tion rate this is clearly the case for most of the domain. We 
can therefore conclude that dynamo action should be possible 
in that case. However, a more detailed comparison would re- 
quire using test fields with wavenumbers other than k - k\, 
as was done in Paper II for a case with slower rotation than 
here. We should therefore not over-interpret our present com- 
parisons. 

The significance of an overshoot layer in connection with 
large-scale dynamo action is not yet entirely clear Tobias 
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Fig. 14. — Vertical profiles of for the same runs as in Fig. [13] The 
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et al. (2008) argue that its presence is unimportant for pro- 
moting the generation of large-scale fields, while Browning 
et al. (2006) find it to be crucial. However, as discussed in 
the present paper, the rotation rate used in Tobias et al. (2008) 
may be too slow for a successful large-scale dynamo. Another 
issue might be the absence of shear-induced helicity fluxes 
(Vishniac & Cho 2001), as discussed in Paper 1. In Browning 
et al. (2006), on the other hand, the presence of shear in the 
overshoot layer plays an important role, as is demonstrated by 
the striking difference to earlier results of Brun et al. (2004) 
without overshoot. Another aspect of the problem is that in lo- 
cal convection simulations the amount of overshoot decreases 
as rotation increases (Ziegler & Riidiger 2003; Kapyla et al. 
2004). 

3.4. Mean-field models 

In order to check how well the turbulent transport coeffi- 
cients presented in the previous section describe the excita- 
tion of large-scale dynamos in the direct simulations, we use a 
one-dimensional mean-field dynamo model in which the test- 
field results for aij and rjij are used as input parameters. The 
model solves the equation 

Ai - aijBj - (rjij + T]6ij)iuoJj, (24) 

where the dot refers to a time derivative. The mean magnetic 
field and current density are given hy B - {-Ay,A^,Q) and 

^oJ - , respectively. Here primes denote z-derivatives. 

The test field results for a/j and rjij can now be directly used 
in the mean-field model, leaving little freedom. We must, 
however, bear in mind the limitations of this simple model. 
Firstly, the transport coefficients were determined for a single 
value of k, whereas in the direct simulations a larger set of 
wavenumbers are available. Secondly, the model is restricted 
to fields that depend only on z, whereas the direct simulations 
indicate that the large-scale field can also vary in the horizon- 
tal directions. 

Nevertheless, it is interesting to study whether the transport 
coefficients derived using the test field procedure can excite a 
dynamo in the present case. The results for the growth rate 
A of the magnetic field from the mean-field models, using the 
coefficients shown in Fig. [13] are shown in Fig. [15] We find 
that the two most rapidly rotating runs show clear dynamo 
action which is consistent with the direct simulations. These 
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Fig. 1 5. — Growth rate of the large-scale magnetic field for the same runs as 
in Fig. ll3l from the one-dimensional dynamo model. 



results support the notion that the mechanism generating the 
large-scale magnetic fields in the direct simulations is the tur- 
bulent (T-effect. 

It should be noted that for values of A different from zero 
the actual growth rate will not be exactly equal to A, because 
in general the components of aij and rjij depend on frequency 
(and thereby on A) - just as they also depend on the value of 
k; see Hubbard & Brandenburg (2009) for details. 

4. CONCLUSIONS 

We use numerical simulations to demonstrate that rigidly 
rotating convection can generate a large-scale dynamo. The 
importance of this result lies in the fact that, according to 
mean-field theory, rotating turbulent convection leads to an 
c-effect, which should result in a large-scale dynamo. So far 
this has not been seen in simulations without shear, leading 
to doubts of the applicability of the mean-field approach. The 
present results show that the lack of large-scale dynamos in 
rotating convection previously reported is probably due to too 
slow rotation. We find that an appreciable large-scale field is 
generated in the direct simulations if Co > 4. 

With the smallest system size explored here, the scale sepa- 
ration is small even in the most rapidly rotating case (Run A6), 
i-C- ^max/^i ~ 2. In this case the large-scale magnetic field is 
concentrated in the k = mode. Doubling the domain size 
increases the scale separation to ^max/^i ~ 5 (Run B6). Here 
the k = contribution is smaller and most of the magnetic 
energy is found for k/ki = 1 and 2 modes. Doubling the size 
of the box once more, k^^^/ki x 10. The magnetic energy 
peaks at the largest scales, although the k = mode is sig- 
nificantly weaker than in the cases with smaller system sizes. 
This suggests that for the present parameters there is a max- 
imum size for the large-scale magnetic field structures which 
is larger than 2d and smaller than 4d. 

We compute the turbulent transport coefficients using the 
test field method and find that the magnitude of the a-effect re- 
mains relatively unaffected when rotation increases. However, 
at the same time the turbulent diffusion is severely quenched 
by the rotation. When these coefficients are used in a one- 
dimensional mean-field dynamo model, the two most rapidly 
rotating runs exhibit a growing dynamo in accordance with 
the direct simulations. Although the mean-field model is too 
simple to fully describe the field in the direct simulations, 
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our results seem to validate the test field method further and 
lend support to our interpretation that the large-scale magnetic 
fields observed in the simulations are due to a turbulent a- 
effect. 

On a more general note, the present results also represent a 
nice demonstration of the usefulness of mean-field theory as 
a predictive tool. In fact, as explained in § 13.21 our work was 
motivated by the earlier findings of Paper 11 that the a-effect 
increases and the turbulent difFusivity decreases as the rotation 
rate is increased. This did already suggest the existence of an 
a^-dynamo for sufficiently rapid rotation - a suggestion that 
we have now been able to confirm in this paper. 

It should be noted that the magnetic field structure of mean- 
field dynamos depends crucially on the geometry of the do- 
main and the nature of the boundary conditions. Therefore, 
our results are not directly relevant to astrophysical bodies, 
because their geometry is not Cartesian nor the boundaries 
periodic. Our models can also not be thought of as local re- 
presentations of a star, although it is possible to get some idea 
about the dependence of the transport coefficients on latitude 
(Ossendrijver et al. 2002; Kapyla et al. 2006; Paper 11). Fur- 
thermore, the values of the magnetic Reynolds number are 
obviously not in the astrophysically relevant regime. Nev- 



ertheless, our simulations provide the first concrete evidence 
of dynamo action from rotating convection without the addi- 
tional help of shear This is an important point, because it 
proves for the first time that an dynamo from rotating con- 
vection exists and that it is strong enough to produce large- 
scale fields. This was thought impossible until now (Cattaneo 
& Hughes 2006; Hughes & Cattaneo 2008; Hughes & Proctor 
2009). 
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